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We study how energy and heat are transferred between two two-dimensional carrier systems 
due to interlayer carrier interaction. We adopt perturbation theory and fluctuation-dissipation 
relations and derive general expression for the interlayer energy transfer rate, which is applicable 
to semiconductor and graphene based bilayers. Our detailed calculations for GaAs and Si bilayers 
suggest that the interlayer energy transfer can exceed the electron-phonon energy transfer below 
(system dependent) finite crossover temperature. For example, for low mobility Si bilayer with 200 
nm layer separation we find a crossover temperature of ~ 1.4 K. 
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Solid bodies that are separated by a small vacuum gap 
can exchange energy and momentum by various mech- 
anisms. When the different bodies represent thermal 
baths we deal with near-field heat transfer and con- 
siderable efforts have been devoted to understand the 
heat transfer via different channels.^ One of the most 
significant exchange channel is built from inter-body 
photon coupling. Surface excitations involving optical 
phonons and plasnions can also play an important role 
and these so-called polariton effects can strongly enhance 
the heat transfer. Recently, also a near-field heat transfer 
channel arising directly from lattice vibrations has been 
proposed.lSl 

Many near-field experiments are in fact performed in 
a configuration where there is no vacuum gap in between 
the different solids. For example, extensive investiga- 
tions of momentum transfer (or drag) due to Coulomb in- 
teraction between two two-dimensional (2D) carrier sys- 





FIG. 1. (color online) (a) Illustration of near-field scatter- 
ing processes (momentum and energy transfer) between 2D 
carrier layers 1 and 2 separated by distance d, at local tem- 
peratures Ti and T2, respectively. A carrier in layer 1 (2) 
experiences momentum scattering fci — > fei + q (fc2 — > fc2 — q) 
due to interlayer interaction. During the process energy AE 
is transferred between the layers, (b) The equivalent ther- 
mal circuit. Pl is the input heating/cooling power (L = 1, 2) 
and G12 is the interlayer thermal conductance. Due to the 
interlayer scattering processes G12 7^ and power (or heat) 
P12 flows between the layers. Layers couple to phonon bath, 
which is at temperature Tp, via electron-phonon thermal con- 
ductance Glp and power Plp flows to the bath. 



tems, that reside in the same heterostructure, have been 
performed.*^ In this work, we investigate how energy and 
heat are transferred between the layers of such bilayer 
carrier system due to interlayer Coulomb interaction (see 
Fig. [ij . In contrast to the vacuum gap case now the car- 
riers couple to the same phonon bath and this coupling 
competes with the interlayer energy transfer [see Fig. 
[ijb)]. We derive general expression of interlayer energy 
transfer rate, which is applicable to semiconductor and 
graphene based bilayers, using perturbation theory and 
fluctuation-dissipation relations. Our formula for inter- 
layer thermal conductance, G12, shows strong similarity 
with the drag resistance formulas that are derived with-in 
similar theoretical framework. Approximation formulas 
and detailed calculations of G12 for semiconductor based 
bilayers are presented. It is shown that 6*12 can dominate 
over electron-phonon thermal conductance. 

The scattering events depicted in Fig. [ija) are medi- 
ated by interlayer interaction which is described by ma- 
trix element Mq. The interlayer Hamiltonian H is given 
by 
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where A is the area, is the electron density operator 
for layer L = 1,2 and c"^\^ is the electron annihilation 
(creation) operator. Variables fc, s and a are wavevector, 
band index and spin index, respectively (here we will 
assume spin degeneracy). All electron variables depend 
on the layer index L, but this is typically not written 
explicitly (e.g. fc = feL)- Factor F^^s is defined by the 
wavefunction of the single particle states and product 
^k-q s'^k,s defines a band form factor. For an ideal 2D 
electron system we have Fk,sa- = 1 and summation over 
band indices s, s' can be ignored. For graphene we have 

-1 denote 



Fk,s = 75 ( 1 se*^*" )^ where s = +1 and s 
conduction and valence bands, respectively, and 9k = 
arctan (ky/kx). 

Next H will be considered as a perturbation Hamilto- 
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nian that will cause transitions from initial state |zl, i2) = 
\i2) with energy Ei = En + i?2i to final state 
1/1, /2) = 1/1) 1/2) with energy Ef = E^j + E^j. Here 
\iL) {\fL)) is the initial (final) state of layer L. The tran- 
sition rate T fi from initial state i to final state / is given 
by the golden rule formula 

T/. = ^ I (/2, /1| H |zl, t2)f 5{E., ~Ef). (2) 

By multiplying T f i by the energy change A_Ei = En — 
El f and performing an ensemble average over the initial 
electronic states, and summing over final electronic states 
we obtain energy transfer rate (heat transfer) 
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where WLi is the weighting factor of carrier layer L 
in state i. We assume that each layer L can be 
described by a local temperature and, therefore, 
WLf = WLiexp[(£;j - Ef)/kBTL]. By using identity 
S{Ea + Eb) = hj^'^duj5{EA~huj)5{EB+huj) and 
definition of correlator 

Chiq-,^) = 2Trh^WLn \ {nL \p^g\ mL)\^ 

n.m 

x5{ELn- EL,n+huj) (4) 

we find 
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Adopting the fluctuation dissipation relatiorP (1 — 
g-Sa;/fcBTi,-)(^^^q^^-) ^ -2hAlm{xL{q,i^)} and the prop- 
erty C[,{q, —oj) = e~''"/'^^'^^Ci(q', w) we find the general 
expression for the interlayer heat transfer 
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where n]^{huj) — {exp{hLu/kBTL) — 1) and Xl{Qj'^) 
the susceptibility, which can depend on T^. At the limit 
7i , T2 — >■ T it is useful to define the interlayer thermal 
conductance GuiT) = Pi2/{Ti ~T2). From Eq. ^ we 
find 



Gi2(r) = 
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which has a striking similarity with bilayer drag resis- 
tance formulaP Equation Q has only single temperature 



and, therefore, it is more convenient to adopt in the case 
studies below. 

We assume that the interlayer interaction is mediated 
by screened Coulomb interaction, when the matrix el- 
ement is given by Mq = ej^2^(q, aj)J7(q)Fi2((i), where 
U{q) = I2ehq is the 2D Fourier transform of Coulomb 
potential (et, is the background dielectric constant) and 
F\2{d) is the spatial form factor, which depends on the 
spatial extent of the electron wave functions and layer 
distance d. For the sake of simplicity here we assume van- 
ishing extent and as a result we have Fi2{d) = exp(— qd). 
The inter-layer dielectric function £12(9, is given 
by £12(9,0;) = [l-U{q)xiiq,u})][l-U{q)x2iQ,^)] - 
F^2U{q)^Xi{q,oo)x2{q,^)' 

In the ballistic limit the carrier mean free path 1^ ex- 
ceeds the layer distance (/g 3> ) and we use the suscep- 
tibility of an ideal 2D electron ga^ 

Xl(9,^) = i'i2zy^ {2z - n^{z,u) - n+{z,u) 

+r-{z,u)~ r+{z,u)}. (8) 



Here Q±{z,u) ~ C±\/{z±u) —1, r±{z,u) = 



iD±yl — {z±u) , z = q/2kp, u — uj/qvp^ C± — 
[z ± u) / \ z ± u\ and D± = for |z ± u| > 1, and C± — 
and D± = 1 for |z ± u| < 1. Here vp (kp) is the Fermi 
velocity (wave vector) and v — v{sp) is the density of 
states at Fermi level ep ^ /ssT. For two similar layers 
we find asymptotic low-temperature result 
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where a ~ 1.86, n — is the screening wave vector 
and / (a) ~ 13.66a~2. The above Equation is valid when 
kpdkBT < 2Ep. Parameter nd characterizes the screen- 
ing of the interlayer interaction: large (small) nd means 
strong (weak) screening. 

In the diffusive limit (wr, le/d ^ 1) the susceptibility 
can be approximated as 



xdq,^) 
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where D = VpT/2 is the diffusion coefficient and r = 
le/vp is the momentum relaxation time. By substituting 
Eq. ([To]) into Eq. ^ we find analytical approximation 
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Here a — e^vD is the DC conductivity of single layer 
and An — T{n)C,{n) = J dxx"~^ / [exp(x) — 1]. Equation 
([n]) is valid when {l^/vp)kBT/h < {k/df. Note that 
Eq. (11) applies also for low mobility graphene bilayers. 

As depicted in the thermal circuit of Fig. [l]jb) G12 
competes with the electron-phonon thermal conductance 
Gpp, which at the limit T^, Tp — T is given by Glp{T) = 
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PlpI {Tl — Tp). In 2D electron systems at low temper- 
atures the electron-phonon energy transfer is dominated 
by screened deformation potential (G£^) and piezoelec- 
tric (Gf^p) interaction with total thermal conductance 
Glp{T) = G^^{T) + Gl^{T). For the deformation po- 
tential contribution we 



(/„(0)S2)z;5:-6t6-", 
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where n = 0(1) represents ballistic (diffusive) limit of 
electron-phonon coupling, for which we have qxrle > 1(< 
1). Here qyr = ksT/fivx is the thermal phonon wave vec- 



vl{t) is the longitudinal 



tor, factor Fi — ttt. ft 

(transversal) phonon velocity and p is the mass density of 
the crystal. The brackets (• • • ) stand for solid angle av- 
erage and 9 is the angle with respect to the z-axis, which 
is perpendicular to the layer plane. (^fn{9)'E?) is an ef- 
fective deformation potential coupling and /o(^) = sin0 
and fi{0) = a+siS e ■ Parameter a — [Klf,VF /vx)~'^ and 
as a result fi{9) sa 1. The piezoelectric coupling gives 
rise to contributiorP^^ 



{U9)K^)vl-^T'-' 
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where (fn{9)K^) is the effective piezo couphng. 

Figure [2] shows Gi2{T) obtained numerically from Eqs. 
([t]) and ([sf in the case of symmetric high mobility (bal- 
listic) GaAs-AlGaAs bilayei^ (here referred as GaAs 
bilayer) with single layer electron density n = 10^^ 
m^2 and d = 20 nm. Asymptotic limit formula of 



In the phonon contribution we 



where 



is also plotted 
.A (sin6'S2)?;-s = iTrS^u^ , wnere 
is the dilatational deformation potential constant 

.-4 e \ „2 _ / 89 ..-4 



10 eV 
and 

^;^(sin(?X2) ^ (^JL^ e2^^(_§_,^ 

where ei4 = —0.16 C/N is the only non-zero element 
of the piezotensor. Oth er pa ram eters can be found from 
Ref.'^. Equations ( 12 ) and ( 13 ) are plotted as symbols 
in Fig. [2] in the ballistic limit. Below few Kelvin piezo- 
electric coupling fully dominates and as a result the tem- 
perature regime where Glj} < Gi2is pushed towards rela- 
tively low, but experimentally achievable, temperatures. 
The crossover occurs at T ~ 140 mK. 

For silicon based bilayeJl^l we consider both ballistic 
and diffusive limits at electron density n — 5 x 10^^ 
m~2. Parameters for Si can be found from Ref.mJ. The 
curves in Fig. [3] are calculated for symmetric high (low) 
mobility Si bilayer system with mobility fj, = 2.5 (0.2) 
m^/Vs, mean-free path 1^ = 200 (16) nm and layer 
distance d = 20 (200) nm. For the high mobility device 
we have used the ballistic limit response [Eq. and 
for the low mobility one the diffusive response [Eq. 
( 10 )] . Even though the layer distance is one magnitude 
larger in the diffusive case, G12 is still significantly larger 
than in the ballistic case. Thus, disorder can strongly 
enhance the interlayer energy and heat transfer. Silicon 
is not piezoelectric so for Glp we need to consider only 
G^p . Due finite electron mean free path (even for 
the high mobility device) we will include ballistic and 
diffusive limits of G^^^f . For simplicity we plot G£^ 
so that it changes abruptly from diffusive to ballistic 
formula [note that Eq. (12 1 is valid close to qxxle — !]■ 




FIG. 2. (color online) The interlayer thermal conductance 
(G12) and deformation potential {Glp) and piezoelectric 
(G^p) electron-phonon thermal conductance at the ballistic 
limit for GaAs bilayer system with electron density n = 10^^ 
m~^. Solid curve is the result from numerical integration us- 
ing the ballistic response function. Dashed curve obtained 
using Eq. ([9|. 




FIG. 3. (color online) The interlayer thermal conductance 
(G12) and deformation potential (Glp) electron-phonon ther- 
mal conductance for two Si bilayer devices with n = 5 x 10^^ 
m^^. For device A (B) solid curves is a result from numeri- 
cal integration using the ballistic (diffusive) susceptibility and 
dashed curves are obtained from the approximation formula 
of Eq. (|9| (Eq. ([lT|). The circle marks the crossover where 



electron-phonon coupling changes from ballistic to diffusive. 
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where S„ is the uniaxial deformation potential constant. 
We use the typical values Sj;(„) = —11.7(9.0) eV. For 
the high and low mobility Si systems the crossover 
temperature where G12 = Glp is 660 mK and 1.4 K, 
respectively. 



Another widely explored bilayer carrier system, which 
can be realized using compound semiconductors^- or 
Sp^, is the electron-hole bilayer. The complexity of the 
valence band makes this system more difficult to analyze 
theoretically. We will not present G12 for such system 
here explicitly, but it should behave in similar fashion 
as its electron- electron counter part. However, one thing 
that may differ drastically from electron-electron bilayer 
system is the carrier-phonon coupling. Due to asymme- 
try in the deformation potential coupling between the 
different layers the carrier-phonon coupling can be un- 
screened and as a result Glp can be strongly enhanced 
at low temperatures.* The enhancement factor depends 
on the details, but in many cases it is of the order of 
{k/qxt)'^, which suggests that for electron-hole bilayers 
Glp can dominate over G12 even down to very low tem- 
peratures. Note that also in symmetric electron bilayers 
Glp can be affected by the presence of another carrier 
system in non-trivial way, but significant enhancement is 
not expectecP. 

The interlayer heat transfer can be investigated experi- 
mentally by varying the input powers Pl while measuring 
the electron temperatures Tl [see Fig. ^h)]. Uniform 
input power follows, for example, from Ohmic heating. 
The electron-phonon contributions G^p can be investi- 
gated independently from G12 by depleting the other car- 
rier system and/or at balanced input power that gives 
Ti = T2. It is not necessarily trivial to measure elec- 
tron temperature of the individual layers. For exam- 



ple, quantum corrections of resistivity and Shubnikov- de 
Haas oscillations, that have been used as electron ther- 
mometer, can be affected by the other layer in non-trivial 
way. More local temperature probes based, e.g., on quan- 
tum point contacts and quantum dots^'^ have also been 
demonstrated. Especially in the case of Si, one possibil- 
ity is to control P/, and measure temperature in doped 
contact regions, when G12 can be investigated by ad- 
justing layer densities by gate voltage. Electron-phonon 
coupling in doped semiconductors can be relative weak 
so that G12 dominates. However, it is important to note 
(as already pointed out in Ref.^l) that similar near-field 
thermal coupling to G12 can exist between bilayer car- 
riers and external gate electrodes. The sign of Pl can 
be also reversed which is equivalent to cooling. This 
can be achieved by quantum dots^^ or semiconductor- 
superconductor contactJi^. In general, G12 could be used 
as a controlled thermalization path (thermal switch) at 
low temperatures. 

In summary, by using first order perturbation theory 
and fluctuation-dissipation relations we derived a gen- 
eral expression of inter-layer energy transfer rate [Eq. 
^] and thermal conductance [Eq. ([t])]. Our formula- 
tion can be applied, for example, to semiconductor and 
graphene based bilayers. We presented analytical ap- 
proximation formulas and detailed calculations for GaAs 
and Si based bilayers. Interlayer energy transfer was com- 
pared to electron-phonon coupling and our results showed 
that the inter-layer energy transfer exceeds the electron- 
phonon energy transfer below crossover temperature of 
140 mK (660 mK) for ballistic GaAs (Si) bilayer with 
d = 20 nm layer distance. Interlayer energy transfer is 
enhanced in the diffusive limit and for low mobility Si 
bilayer with d — 200 nm the crossover occurs at ^ 1.4 K. 
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